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Abstract 


The  modeling  of  a  moving  target  for  a  single  channel  spaceborne  SAR  geometry  has  al¬ 
ready  been  accomplished  to  a  high  degree  of  accuracy  by  Eldhuset  et  al. ,  but  extending  the 
model  to  include  a  SAR  system  that  is  equipped  with  multiple  apertures  (e.g.  RADARSAT- 
2,  TerraSAR-X)  still  requires  further  work.  The  purpose  of  this  short  technical  memo  is 
to  do  exactly  that  -  to  derive  a  set  of  equations  of  motion  of  a  ground  moving  target  for  a 
multi-channel  spaceborne  SAR.  These  equations  of  motion  will  be  shown  to  be  applicable 
to  both  airborne  and  spaceborne  multi-channel  SAR  systems  in  stripmap  mode. 


Resume 


Eldhuset  et  coll,  ont  deja  modelise  avec  grande  exactitude  la  geometrie  de  l’observation 
d’une  cible  mobile  par  un  RSO  monocanal  spatioporte.  Toutefois,  de  nouveaux  travaux  sont 
necessaires  pour  adapter  leur  modele  a  un  RSO  a  ouvertures  multiples  (p.  ex.,  RADARS AT- 
2  ou  TerraSAR-X).  Les  travaux  menant  au  present  memoire  technique  visaient  a  deriver  un 
ensemble  d’ equations  de  mouvement  d’une  cible  terrestre  mobile  suivie  par  un  RSO  multi- 
canal  spatioporte.  Nous  montrons  que  Ton  peut  appliquer  ces  equations  a  l’observation  en 
bande,  autant  par  un  RSO  multicanal  aeroporte  que  spatioporte. 
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Executive  summary 


Equations  of  motion  of  a  ground  moving  target  for  a 
multi-channel  spaceborne  SAR 

S.  Chiu,  M.V.  Dragosevic;  DRDC  Ottawa  TM  2008-326;  Defence  R&D  Canada  - 

Ottawa;  March  2009. 

Background:  Both  high  resolution  Synthetic  Aperture  Radar  (SAR)  processing  and  SAR 
Ground  Moving  Target  Indication  (SAR-GMTI)  require  that  a  highly  accurate  imaging 
geometry  model  be  first  established.  This  can  be  quite  easily  accomplished  for  the  case  of 
an  airborne  platform,  which  is  assumed  to  be  moving  along  a  straight  line  and  transmitting 
uniformly  spaced  pulses.  This  assumption  requires  good  motion  compensation  and  good 
control  of  the  pulse  repetition  frequency  (PRF)  as  a  function  of  ground  speed.  For  an 
orbiting  platform,  the  gravitational  force  plays  a  key  role  in  defining  the  satellite  trajectory 
and,  thus,  must  be  taken  into  account.  Also,  equations  of  motion  of  a  moving  target  that 
accurately  model  a  SAR  system  equipped  with  multiple  apertures  (e.g.  RADARSAT-2, 
TerraS  AR-X)  are  evidently  absent  in  the  open  literature,  partly  because  there  were  no  multi¬ 
aperture  spaceborne  SARs  in  the  unclassified  world  prior  to  late  2007. 


Principal  results:  The  purpose  of  this  short  technical  memo  is  to  derive  a  set  of  equations 
of  motion  that  accurately  describes  a  ground  moving  target  in  spaceborne  multi-channel 
SAR  imaging  geometry.  A  suitable  multi-channel  model  of  the  target  range  history  is 
derived  in  the  Earth  Centered,  Earth-Fixed  (ECEF)  frame  of  reference  using  linearization 
for  small  angles  as  a  function  of  the  receive  phase  center  index  p  (=  1,  2,  3,  etc.)  along  the 
antenna  and  the  pulse  (or  slow)  time  t: 


RJt)  =R0  +  (p-l)D('¥- d>0)  +  (Vtr  -  Fsr)  (t  - 10) 


+ 


(P~1)D 

Ro 


UVl~2W.+Atr)  f 

Ro 


The  above  range  equation  includes  the  following  system  parameters:  Ro  is  the  target  range 
at  arbitrary  time  to;  D  is  the  distance  between  the  effective  phase  centers;  'F  is  the  so- 
called  squint  angle  between  the  radar  broadside  in  the  slant-range  direction  and  the  antenna 
pointing  vector;  d>o  is  the  direction  of  arrival  angle  (measured  off-broadside  from  the  SAR); 
Vtr  and  14,  are  the  radial  velocity  components  of  the  target  and  the  SAR,  respectively;  Vta  is 
the  target  velocity  component  in  the  along-track  or  flight  direction;  14  and  Vg  are  the  SAR 
orbital  and  ground  speeds,  respectively;  14  is  the  so-called  effective  velocity  of  the  satellite; 
and  Atr  is  the  target  radial  acceleration  component.  Ro.  d>o,  Vtr,  Vta,  and  Atr  are  defined  for 
time  to.  14,  depends  on  14  and  <t>o  in  a  predictable  way.  14,  Vg,  and  14  vary  slowly  with  time 
and,  therefore,  may  be  evaluated  anywhere  in  the  neighborhood  of  to. 
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Significance  of  results:  The  accuracy  of  the  equations  of  motion  derived  in  this  short 
memo  has  been  tested  and  validated  using  the  recently  acquired  RADARSAT-2  MODEX 
data.  These  equations  of  motion  have  been  shown  to  be  applicable  to  both  airborne  and 
spaceborne  stripmap  imaging  geometries  and  serve  as  a  physical  basis  for  further  algorithm 
development. 
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Equations  of  motion  of  a  ground  moving  target  for  a 
multi-channel  spaceborne  SAR 

S.  Chiu,  M.V.  Dragosevic ;  DRDC  Ottawa  TM  2008-326  ;  R  &  D  pour  la  defense 

Canada  -  Ottawa ;  mars  2009. 

Introduction  :  Le  traitement  des  donnees  des  radars  a  synthese  d'ouverture  (RSO)  de 
haute  resolution  et  1’ indication  de  cibles  terrestres  mobiles  au  sol  necessitent  Petablisse- 
ment  prealable  d’un  modele  tres  precis  de  la  geometrie  de  la  visee.  Ce  modele  est  facile- 
ment  realise  pour  un  RSO  aeroporte  pour  lequel  on  peut  postuler  une  trajectoire  rectiligne 
et  remission  reguliere  d’ impulsions.  Pour  etre  valides,  ces  deux  hypotheses  necessitent  une 
bonne  compensation  de  mouvement  et  une  regulation  adequate  de  la  frequence  d’ emission 
des  impulsions  en  fonction  de  la  vitesse  au  sol.  Dans  le  cas  d’une  plateforme  orbitale,  la 
gravitation  est  un  element  essentiel  de  la  definition  de  la  trajectoire  et  nous  devons  en  tenir 
compte.  Or,  dans  les  sources  publiees,  on  ne  retrouve  pas  d’ equations  de  mouvement  d'une 
cible  mobile  permettant  de  modeliser  avec  exactitude  un  systeme  RSO  a  ouvertures  mul¬ 
tiples  (p.  ex.,  RADARSAT-2  ou  TerraS AR-X),  notamment  parce  qu’avant  la  fin  de  2007,  il 
n’existait  aucun  systeme  satellise  semblable  dans  le  domaine  non  classifie. 


Resultats  :  Le  present  memoire  technique  presente  notre  derivation  d’un  ensemble  d’ equa¬ 
tions  de  mouvement  decrivant  avec  exactitude  une  cible  terrestre  mobile  observee  par  un 
RSO  multicanal  spatioporte.  Pour  un  repere  geocentrique  a  axes  fixes,  nous  avons  derive 
un  modele  adequat  de  P  evolution  des  distances  entre  la  cible  et  un  RSO  multicanal,  par 
linearisation  des  petits  angles,  en  fonction  de  l’indice  p  (=  1,  2,  3,  etc.)  du  centre  de  phase 
de  reception  le  long  de  l’antenne  et  le  temps  t  de  P impulsion  (lent)  : 


Rp(t)  =  R0  +  (p  - \)D  OP  -  d>0)  +  (Vtr  -  Vsr)  ( t  - 10) 


+ 


(p-1)£> 

Ro 


<yu-vs)(t-to)+U 


'^-2V,Vt»+AiA(t_io)2 


R 


0 


Les  equations  ci-dessus  pour  la  distance  comportent  les  parametres  suivants  :  Rq  est  la 
distance  de  la  cible  au  temps  arbitraire  to-  D  est  la  distance  entre  les  centres  effectifs  de 
phase.  'P  correspond  a  l’angle  de  concentration  entre  la  position  de  rayonnement  transverse 
du  radar  dans  le  plan  inclinaison-temps  et  le  vecteur  de  pointage  de  Pantenne.  4>o  represente 
la  direction  de  P  angle  d’ incidence  de  l’echo  (mesure  depuis  le  cote  du  RSO).  Vtr  et  Vsr 
sont  respectivement  les  vitesses  radiales  de  la  cible  et  du  RSO.  Vta  est  la  composante  de 
la  vitesse  de  la  cible  le  long  de  la  bande  observee  ou  de  la  trajectoire  de  vol.  Vs  et  Vg  sont 
respectivement  les  vitesses  du  RSO  en  orbite  et  au  sol.  V e  est  la  vitesse  effective  du  satellite. 
Atr  est  la  composante  radiale  de  P acceleration  de  la  cible.  Les  valeurs  Rq,  d>o,  Vtr,  Vta,  et 
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v 


Atr  sont  definies  pour  le  temps  to-  V sr  depend  de  V s  et  <5>o  de  I’acon  previsible.  Puisque  Vs, 
Vg,  et  V e  varient  lentement  en  fonction  du  temps,  on  peut  les  evaluer  aux  environs  de  to. 

Portee  :  Nous  avons  teste  et  valide  1’ exactitude  des  equations  de  mouvement  qui  figurent 
dans  le  present  memoire  a  l'aide  de  donnees  RADARSAT-2  recemment  acquises  lors  des 
essais  MODEX.  Nous  avons  demontre  qu’elles  peuvent  s’appliquer  a  la  geometrie  de  l’ima- 
gerie  en  bande  (mode  stripmap)  par  un  radar  aeroporte  ou  spatioporte.  Elies  constituent 
egalement  un  fondement  tangible  pour  la  mise  au  point  d’autres  algorithmes. 


VI 
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1  Introduction 


High  resolution  Synthetic  Aperture  Radar  (SAR)  processing  requires  that  a  highly  accurate 
imaging  geometry  model  be  first  established.  For  SAR  Ground  Moving  Target  Indica¬ 
tion  (SAR-GMTI),  the  underlying  assumption  that  the  radar  scene  is  stationary  must  be 
extended  to  include  non-stationary  scenes  or  moving  targets.  This  can  be  quite  easily  ac¬ 
complished  for  the  case  of  an  airborne  platform  [1],  which  is  assumed  to  be  moving  along  a 
straight  line  and  transmitting  uniformly  spaced  pulses.  This  assumption  requires  good  mo¬ 
tion  compensation  and  good  control  of  the  pulse  repetition  frequency  (PRF)  as  a  function 
of  ground  speed.  The  same  cannot  be  said  about  a  spaceborne  platform,  where  the  earth’s 
gravitational  force  plays  a  key  role  in  defining  the  platform  trajectory  and  the  velocity  of 
the  radar  antenna  footprint  as  it  sweeps  along  the  surface  of  the  earth.  The  modeling  of 
a  moving  target  for  a  single  channel  spaceborne  SAR  geometry  has  already  been  accom¬ 
plished  to  a  high  degree  of  accuracy  by  Eldhuset  et  al.  [2]  [3].  However,  the  extension  of 
the  model  to  include  a  SAR  system  that  is  equipped  with  multiple  apertures  is  evidently  ab¬ 
sent  in  the  open  literature,  partly  because  there  were  no  existing  spaceborne  SAR  systems 
in  the  unclassified  world  equipped  with  such  a  capability  until  the  launches  of  TerraS AR-X 
and  RADARSAT-2  in  late  2007.  The  purpose  of  this  short  technical  memo  is  to  derive  a 
set  of  equations  of  motion  of  a  ground  moving  target  for  a  multi-channel  spaceborne  SAR. 
These  equations  of  motion  will  be  shown  to  be  applicable  to  both  airborne  and  spaceborne 
stripmap  imaging  geometries  and  serve  as  a  physical  basis  for  further  algorithm  develop¬ 
ment. 


2  Equations  of  Motion  of  a  Moving  Target 


The  relative  position  vector  of  a  moving  target  with  respect  to  an  imaging  SAR  satellite,  in 
the  earth  centered  earth  fixed  (ECEF)  system,  can  be  written  as 

R,  =  Rf  —  R^ ,  (1) 

where  indices  t  and  s  denote  ‘target’  and  ‘satellite,’  respectively,  and  bold  letters  indicate 
vectors.  The  magnitude  of  a  vector  is  represented  by  the  same  symbol  with  the  regular  italic 
font.  In  the  ECEF  frame,  the  earth  motion  is  absorbed  into  the  relative  satellite  motion. 

The  Doppler  centroid  and  Doppler  rate  are  proportional  to  R  and  R,  respectively,  where  the 
dot  and  double-dot  notations  indicate  first  and  second  derivatives  with  respect  to  time.  A 
common  approach  to  the  derivation  of  R  and  R  is  to  start  from  the  identity  [3] 

R2  =  RrR  ,  (2) 

and  to  differentiate  it  with  respect  to  time.  (•)r  denotes  the  vector  (or  matrix)  transpose. 


DRDC  Ottawa  TM  2008-326 


1 


Differentiating  both  sides  of  (2)  with  respect  to  time,  we  get 

(3a) 
(3b) 

where  Ay  B  =  BrA  is  obviously  true  for  any  two  vectors  A  and  B. 


2  RR  =  RrR  +  RrR 
RrR  RrR 

R  = 


R 


R 


Equation  (3b)  can  be  rewritten  as 

Rr 

R  =  Rs) 

-  -^(Vr-v,) 

=  Vtr-Vsn 


(4a) 

(4b) 

(4c) 


where  Vr  =  Rf  is  the  velocity  vector  of  the  moving  target,  Vs  =  Rs  is  the  velocity  vector 
of  the  satellite,  and 


Vtr  = 
vsr  = 


(5a) 

(5b) 


are  the  projections  of  the  target  and  satellite  velocity  vectors  onto  the  line  of  sight  (LOS) 
or  the  radial  direction,  respectively.  Also,  the  radial  speed  of  a  stationary  target  as  “seen” 
by  the  radar  due  to  the  platform  motion  is  equal  to  —Vsr.  Therefore,  the  Doppler  shift  at 
the  beam  center  induced  by  the  motion  of  the  platform  (or  the  stationary  clutter  Doppler 
centroid)  is  given  by 

Vsr 

/dc  —  2-£-  ,  (6) 

and  the  Doppler  shift  due  to  the  target’s  radial  speed  is 

Vtr 

fdc  =  —2-35  .  (7) 

Therefore,  the  total  Doppler  shift  is  given  by 


Fdc  =  —2 


Vr 

X 


(8) 


Again,  differentiating  both  sides  of  (3  a)  with  respect  to  time  yields 


2  (R2+RR) 

=  RrR  +  RrR  +  RrR  +  RrR 

(9a) 

=  2RrR  +  2RrR 

(9b) 

RR 

=  RrR  +  RrR -R2  . 

(9c) 

2 
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Using  the  following  definitions 


V, 

=  R , 

(10a) 

V, 

=  R 

(10b) 

Af 

=  Rr 

(10c) 

A, 

=  Rs 

(10d) 

A 

Oi 

< 

1 

<>> 

1 

III 

(10e) 

Atr 

=  RrA  t/R, 

(10f) 

(9c)  can  be  rewritten  as 

RR  =  -RrA,  +  (Vr-V,)r(Vf-V,)-(U,--UT)2  (Ha) 

=  RAtr -  RrAs  +  V2  -  Vf V,.  -  V[ V,  +  V2  -  (V?r -  2V,rVsr  +  V2 )  (lib) 
=  V2  -  Rr A,  +  RAtr  +  Vf  -  2VS  (£V,)  -  Vfr  +  2VtrVsr  -  V2  .  (11c) 


Therefore, 


where 


R 


Vl_Vl+A  |  V?  2VsVta  jg  2V^. 

7?  R  fr  R  7?  R  R 


(12) 


l/2  =  V2-RrA,  (13a) 

/VAr 

R,  =  -/  V,.  (13b) 

\  RS  / 


Ve  is  the  so-called  “effective  velocity”  often  used  in  the  spaceborne  SAR  processing  to 
model  the  range  equation  and  Vta  is  the  projection  of  the  target  velocity  onto  the  direction  of 
platform  velocity  Vs,  which  is  also  called  the  along-track  direction  and  may  not  necessarily 
be  parallel  to  the  ground  track. 


The  instantaneous  slant  range  equation  (or  history)  R{t)  is  the  key  to  high  precision  SAR 
processing.  The  accurate  estimation  of  the  effective  velocity  in  the  equation  allows  compli¬ 
cated  mathematical  manipulations  involving  a  satellite/earth  geometry  model  to  be  avoided 
and  a  simple  hyperbolic  approximation  to  be  adopted  in  most  high  precision  SAR  process¬ 
ing  algorithms.  The  hyperbolic  model  can  be  further  simplified  and  approximated  using  a 
second  order  Taylor  series  expansion  or  a  parabolic  model  without  significantly  incurring 
further  loss  of  accuracy,  considering  typical  RADARSAT-2  dwell  times  and  resolution. 
This  may  not  be  true  in  general. 

If  R  and  R  in  (4a)  and  (12)  are  evaluated  at  some  arbitrary  time  to,  then  the  range  equation 
can  be  approximated  by  the  Taylor  series  expansion: 

R(t)  ~Ro  +  Vro(t  -to)  +  ^(t-  to)2  ,  (14) 
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where 


Rq  =  R(to)  =  |Rf(fo)  —  R-s(^o)| 

Vro  =  R(to)  =  ^-[Vf(fo)  -  Vs(?0)]  =  Vtr -  VS} 
Ro 


A  /-(]  —  R(to)  — 


v}  -  vl  V ;2  -  V2  +  2 VtrVSr  -  2VsVta 


r  sr  _|_  v  t 


V2  -V2 

v  e  y  sr 


Ro 


+ 


Ro 

2VtrVsr-2VsVta 


R 


o 


Ro 

+Atr . 


+  A, 


tr 


(15a) 

(15b) 

(15c) 

(15d) 


Vf  and  V2  are  considered  negligible  with  respect  to  ‘ 2VsVta ’  and,  therefore,  are  dropped  in 
(15d).  If  to  is  chosen  to  be  the  broadside  time  //,,  then  Vsr  =  0  by  definition  and  the  radial 
direction  (subscripted  by  V)  becomes  exactly  perpendicular  to  the  flight  direction  or  the 
along-track  direction  (subscripted  by  ‘a’).  Under  this  condition,  (15b)  and  (15d)  become 

Vrb  =  R(tb)  =  ^\t(tb)=Vtr  (16a) 

Rb 

.  j±/  \  _  Ve  2UsV)a 

A)-b  ~  R{tb)  —  —  h  Afr ,  (16b) 

Rb 


where  Vtr  and  Atr  are  now  the  target’s  down-range  (or  across-track)  velocity  and  accelera¬ 
tion  components,  respectively,  and  Rb  is  the  broadside  range  of  the  moving  target. 


3  Two-Channel  SAR 


The  use  of  a  parabolic  model  is  useful  in  the  derivation  of  range  equations  for  multi-channel 
SAR  systems.  In  this  section  the  range  equation  for  the  second  aperture  of  a  two-channel 
SAR  is  derived. 


3.1  Local  Frame  of  Reference 


In  order  to  continue  with  our  derivations,  we  first  define  a  local  flight  (LF)  frame  of  refer¬ 
ence  for  the  radar  as  shown  in  Figure  1,  where  d  is  defined  as  the  unit  vector  pointing  down 
from  the  radar’s  center  of  gravity  to  the  center  of  the  earth.  To  define  the  second  axis,  we 
vector  cross-multiply  d  with  the  radar’s  velocity  vector  Vs  to  form  the  right  pointing  unit 
vector  r: 


_  dx  Vs 
r  |dxV,r 

Then  the  third  unit  vector,  which  completes  the  local  reference  frame,  is  given  by 


(17) 


f=rxd. 


(18) 


We  should  point  out  that  Vs  is  not  necessarily  in  the  exact  same  direction  as  f ,  as  illustrated 
in  Figure  1 


4 
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Figure  1 :  Local  reference  frame  of  radar. 
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3.2  Transformation  Matrix 


We  now  derive  the  transformation  matrix  from  the  LF  reference  frame  to  the  ECEF  refer¬ 
ence  frame.  To  begin,  we  express  the  unit  vector  d  in  the  ECEF  frame: 


Then  r  becomes 


where 


d 


R-s 

I R® 

1 

Rs 


R 

R 

R 


SX 

sy 

sz 


d  x  V, 

|d  x  V, 

d  x  V, 


V/ior 

1 

RsVhor 


RszVsy  RsyVsz 
Rsx^sz  Rsz^sx 
RsyVsx  RsxVsy 


(19a) 

(19b) 

(20a) 

(20b) 

(20c) 


Vhor 


V, 


d  x  V, 

'  Vnc  " 

Vsy 

Vsz 


(21a) 

(21b) 


Vhor  is  obviously  the  horizontal  velocity  component  of  the  radar  and  can  be  easily  shown 
to  be 


Vhor  — 


+  V2  i  y2  _y2 
1  v  sy  '  v  sz  vet 


(22) 


where  Vver  is  the  vertical  velocity  component  of  the  radar  platform  and  is  given  by 

T/  _  R-J,r  _  R.sxVsx  T  RsyVsy  +  ^.vzCsz 


(23) 


We  are  now  ready  to  express  the  forward  unit  vector  f  in  the  ECEF  frame  as 


f 


rxd= — d  x  r 


1 

Rs 


Rs 

R 

R 


sy 


sz 


1 

RsVhor 


RszVsy  RsyVsz 
RsxVsz  RszVsx 
RsyVsx  RsxVsy 


1 

RsVhor 


RsVsx  RsxYver 
Rs  Vsy  RsyVver 
RsVsz  RszYver 


(24a) 

(24b) 


(24c) 
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Finally,  the  transformation  matrix  from  the  LF  reference  frame  to  the  ECEF  frame  [4]  is 
simply 

Yf  =  [  f  r  d  ]  (25a) 

j  Rs^sx  RsxYver  RszYsy  Rsy^sz  Rsx^hor 

—  ~d~\/  Rs^sy  Rsy^ver  Rsx^sz  RszYsx  Rsy^hor  •  (25b) 

S  h0T  [  RsVsz-RszVver  Rs}VSx-RsxVsy  -RszVhor  . 

3.3  Antenna  Look  Vector 

Let  the  ideal  look  direction  of  the  antenna  in  the  LF  frame,  with  an  off-nadir  angle  0 
pointing  at  a  zero  Doppler  point  on  the  surface  of  the  earth,  be 

0 

uo  =  sin0  .  (26) 

COS0 

Then  the  actual  antenna  look  vector  (or  pointing  vector)  in  the  local  reference  frame  of  the 
radar  is  given  by 

u  =  r^uo 

— cpvsin0  +  cppcos0 
sin0 

COS0 

where  is  the  yaw-pitch  rotation  matrix,  and  cpv  and  (p/;  are  the  yaw  and  pitch  angles 
about  the  axes  d  and  r,  respectively.  We  are  assuming  that  cpv  and  (p/;  correspond  to  a 
LOS  within  the  beam,  but  not  necessarily  at  its  center.  For  RADARSAT-2,  cpv  and  cp;, 
are  typically  small  (<C  1)  in  the  ECEF  frame  due  to  the  zero-Doppler  beam  steering.  The 
rotation  matrix  r(p  can,  therefore,  be  shown  as 

1  -(pv  Cp  p 

Y(p  =  MrMd  =  %  1 

_  -<PP  0  1 

1  —  CPy  9  p  ' 

«  CPy  1  0  , 

.  “<Pp  0  1 

where 

coscpp  0  sincpp  1  0  cp/7 

M,.  =  0  10  ^010 

—  sincpp  0  coscpp  —  cpp  0  1 

cos  cp^.  —  sincpy  0  1  —  cpv  0 

Mf/  =  sincpv  coscpv  0  rs  cpv  1  0 

0  0  1  J  [  0  0  1  _ 


(29a) 

(29b) 


(28a) 

(28b) 


(27a) 

(27b) 
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The  term,  cpycpp,  is  considered  negligible  [4]  and  is  set  to  zero  in  (27b)  and  (28b).  In  the 
ECEF  frame,  the  antenna  look  vector  u  is  then  given  by 

u  =  r/u  =  r/r(pu0 .  (30) 


We  should  note  that  the  look  vector  u  is  not  necessarily  in  the  direction  of  the  beam  center, 
rather  it  points  to  the  direction  of  the  target  of  interest  within  the  beam  footprint. 


3.4  Displacement  Vector 


Let  D  denote  the  vector  pointing  from  the  effective  phase  center  of  the  aft  sub- aperture  to 
the  effective  phase  center  of  the  fore  sub-aperture  in  the  LF  frame,  D  can  then  be  expressed 
as 


where 


D 

1 

D  =  rv 

0 

ftj  I) 

Yy 

0 

.  — Yp  . 

1 

-Yy 

r<p(cp  =  y)  ~ 

Yy 

1 

-Vp  o 


5 


Vp 

0 

1 


(31) 


(32) 


and  \\fy  and  \\lp  are  the  pitch  and  yaw  angles  (or  the  orientation)  of  the  antenna,  representing 
the  attitude  of  the  spacecraft  in  the  LF  frame  of  reference.  In  the  ECEF  frame,  D  becomes 


D 


T  fD  =  DT 


f 


1 

Yy 

-Yp 


(33) 


3.5  Range  Equations  for  Multiple  Phase  Centers 

A  two-aperture  SAR-MTI  system  is  again  assumed  in  the  following  derivations  with  the 
understanding  that  the  derived  equations  can  be  generalized  to  a  multi- aperture  system. 
Let  Rv  |  and  Rs2  denote  the  position  vectors  of  the  antenna’s  two  effective  (or  two-way) 
phase  centers  in  the  ECEF  frame,  respectively.  The  aft  antenna  phase  center  Rs2  is  then 
displaced  from  the  fore  antenna  phase  center  R4i  by  — D.  For  the  case  of  RADARSAT-2, 
the  displacement  vector  D  is  closely  aligned  with  the  radar’s  velocity  vector  Vs.  Perfect 
alignment  would  be  optimal  because  it  would  allow  the  aft  phase  center  to  pass  through 
the  same  ECEF  position  as  the  fore  phase  center  with  a  time  delay  of  D/Vs,  where  D  is 
the  distance  between  the  two  effective  phase  centers.  This  perfect  alignment  would  also 
mean  that  the  whole  antenna  is  ideally  steered,  generating  a  zero  Doppler  centroid  in  the 
clutter  Doppler  spectrum.  In  the  presence  of  a  non-zero  Doppler  centroid,  there  exists  a 
non-zero  across-track  component  of  D,  which  translates  into  a  small  across-track  baseline. 
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In  the  case  of  a  real  spaceborne  SAR-MTI  system,  such  as  the  RADARSAT-2  MODEX, 
this  small  cross-track  component  is  always  present  and,  therefore,  must  be  compensated 
for  or  taken  into  account  in  the  system  modeling. 

The  slant-range  vector  R2  from  the  aft  antenna  phase  center  to  the  target  can,  therefore,  be 
expressed  as 


R  2  =  Rf-  R.s2 

=  Rf  -  (Ryi  -  D) 

=  Rj  +  D  , 


(34a) 

(34b) 

(34c) 


where  Ri  =  Rf  —  Rvi  and  Rv2  =  Rv  1  —  D.  Then  the  projections  of  these  slant-range 
vectors,  Ri  and  R2,  along  the  look  vector  u  direction  are  given  by 


R\  =  Rfu 

R2  =  R(u  =  (Rf  +  Dr)u  =  R\  +  Dru 
=  R\  +  (Ty-Dj^Tyu  =  R{  +  D^rjryii 
=  R\  +Dru  . 


From  (27b)  and  (31),  (35d)  becomes 


R2(t)  «  /?!(*)+£>[!  XRy  ~Vp' 


— cpvsin0  +  cp/;cos0 
sin0 

COS0 


=  Ri{t)  +R[(\|/y  —  cpy)  sin0  —  (\|/p  —  cpp)  cos0] 
=  Ri(f)+D(^-4>), 


(35a) 

(35b) 

(35c) 

(35d) 


(36a) 

(36b) 

(36c) 


where 


'T  —\\fy  sin  0  —  \| fp  cos  0  (37a) 

=cpv  sin  0  —  cpp  cos  0  .  (37b) 


and  O  are  now  measured  in  the  slant-range  plane.  As  the  antenna  footprint  sweeps 
across  the  target,  the  pitch  angle  cpp  hardly  changes  (i.e.  remains  virtually  constant)  such 
that  cpp  ^  \\ ip,  resulting  in  —  <J>  r;  (v)/y  —  (py)  sin  0.  In  the  case  of  RADARSAT-2,  \) iy  and  \\fp 
are  usually  small  but  non-zero  such  that  the  beam  center  is  not  located  exactly  at  the  zero- 
Doppler  point  on  the  surface  of  the  earth  (in  the  ECEF  frame).  This  residual  beam  squint 
T*  generates  a  small  constant  along-track  interferometric  phase,  which  is  usually  removed 
by  the  digital-balance  processing  of  the  signal  channels  and  can,  therefore,  be  ignored.  For 
the  sake  of  completeness,  however,  we  shall  keep  the  term  in  (36c).  Then,  the  zeroth-order 
coefficient  of  the  Taylor  expansion  of  R2(t)  evaluated  at  arbitrary  time  To  can  be  expressed 
as 

fl2(io)=^o)-£>  [<£>(*)) •  (38) 
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Next,  we  derive  the  first-order  coefficient  of  the  Taylor  series  expansion  of  Riit).  From 
(34c),  we  obtain 


4 

R1R2 


Ri 


(Ri+D)r(Ri+D) 

(Ri+D)r(Ri+D) 

R[Ri  +  R[D  +  DrRi  +  DrD 


R\R\  hi 

+  1 


Ri 

+ 


RfD  DrRi  Drr/D 


+ 


R2  Ri  R2  Ri 

Rr  •  Dr  .  0(D2) 

Ri+itd+~r~  r 

K, +urD+°  I'V.-V.J, 

K 


(39a) 

(39b) 

(39c) 

(39d) 

(39e) 

(39f) 


where  it  can  be  shown  that  R^Ri  =  R\R\,  R\  ~  R2  —  R,  and  the  0(D2)  term  can  be 
neglected. 


First,  we  derive  the  second  term  in  (39f): 

urD  =  nT  jt  (T/D)  =  ur  (tfi>  +  176)  =  urf ft)  ,  (40) 

where  we  have  assumed  that  the  spacecraft  attitude  is  not  changing  in  the  LF  frame  such 
that  time  derivatives  of  \j/v  and  \\fp  (or  D)  are  equal  to  zero  in  the  imaging  time  inter¬ 
val.  We  also  assume,  for  simplicity,  that  l|/v  and  \\rp  are  small.  This  is  normally  true  for 
RADARSAT-2,  which  utilizes  dynamic  antenna  squinting  to  constantly  steer  the  beam  to  a 
zero  Doppler  point  on  the  surface  of  the  earth.  Therefore,  (40)  becomes 


r  1  1 

r  1 1 

urD  =  u  TtfD 

1 

« uTTfD 

- 1 

0  0 

(41) 


Here,  we  need  to  find  the  first  time  derivative  of  Tf  (or  T  /■),  which  can  be  shown  to  be 


1 

RsVhor 


RsVsx  T  RrA  sx  VsxVver  Rsx^ver 

R\ V sy  -j-  RsASy  Vsy^ver  Rsy^ver 
RrVsz  3”  RA  sz.  VszYver  RszYver 


R.sz^sy  R\yA  sz 
Rsx^sz  Rsz^SX 
RsyA  sx  Rsx^sy 


VsxVhor  RsxVho 


VsyVhor  Rsy^hor 
VszYhor  RszYhor 

(42) 


where  terms  of  the  type  VsxVsy,  VSXVSZ,  and  VszVsy  cancel  out  in  the  second  column  of  (42) 
and  are,  therefore,  dropped.  We  can  further  simplify  (42)  by  noting  that  Rs  &  0,  V/wr  &  0, 
and  Vver  ~  0: 


RsV ho 


RrA  sx 

R.A 

RA 


sy  VSyVver 
'  VszYver 


R.SZ.A  sy 
Rs\A  sz 
RsyA  sx  ' 


syAsz 

VsxVhor 

sz.Asx 

—  VsyVhor 

sxA-sy 

VszYhor 

(43) 
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Therefore,  (41)  becomes 


Rs^sx  VsxYver 


LJ  T 

- u 

RsYior 

Rs^sy  VsyVver 

(44a) 

Rs^sz  VszVver 

D  j  . 

T  ,  U  A.y 

DVver  t 

- u  Vs 

(44b) 

Vhor 

RsVhor 

D  j  . 

- u  A, 

Vhor 

(44c) 

The  last  term  in  (44b)  is  very  small,  since  the  look  vector  u  is  virtually  perpendicular  to  Vs 
(i.e.  u  _L  Vs)  in  the  ECEF  frame,  and  is,  therefore,  ignored. 

We  now  derive  the  last  term  of  (39f).  From  (25b)  and  (33), 

Dr 

a-cv.-v.) 


=  p[!  Vv  -¥p]r^  Vt-  Vs 


D[  1  0  o' 

RRs Vhor 

(  \ Vs 

x  Vr  -  V< 


(45a) 


RsVsx  Rsx^ver  RsVsy  RsyVver  Rs^sz  RszYver 

RszYsy  RsyV.sz.  RsxV sz  RszYsx  RsyVsx  RsxV sy 


RsxVhc 


RsyVhor 


RszYhor 


D  RsVsx  RsxYver  RsVsy  RsyVver  R.Ys-  RszYver 

RRsVhor 


(45b) 


(45c) 


where  \|/v  -C  1 ,  \\fp  <C  1,  and  they  are  neglected  in  (45b).  Also  by  noting  that 

RsVs  -Vver  Rj  =  [  RsVsx- RsxV ver  RsVsy  -  RsyVver  RSVSZ  -  RszVver  ]  (46a) 

Rj  =  Rj  T  R;  —  Rf  =  R  T  R/  (46b) 

VjV,  =  VsVta  (46c) 

/rt  i  t">  \T  m  7  i  n  T  a  t  ~  m;  r  mz  x\ 


(R  +  Rf)  Vf  —  RVtr  +  Rr  Vf  ~  RVtr , 


(46d) 
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we  can  rewrite  (45c)  as 

+  DD  — —  [-RsV?  +  Vver  ( RsxVsx  +  RsyVsy  +  ^szVsz)]  (47a) 

KKsVhor 

“  VmR^  v'  “  R’v? + v™Rf  VJ  (47b) 

=  7T7r7r-  kvJV,  -  V„,  (R  +  R,f  V,  -  Rtf  +  V,Ul  (47c) 

KKsVhor  L 

=  -  RVverVtr  -  Rsyj  +  RsV?er)  {Aid) 

RRs^hor 


=k(MiWtMtWJLA  +  (iMM  (47c) 

[\VkorJ  R  \Vh„rJ  R,  \Vhor)  R  \V„or J  R  \  ' 
where  Vta  and  Vtr  are  along-track  and  radial  speeds  of  the  moving  target,  respectively. 

Finally,  (47e)  can  be  further  simplified  by  noting  that  V]wr  ~  Vs  and  Vver/Vflor  <C  1,  yielding 


y(V|  — V,)«^(v„  — r,)  .  (48) 

Putting  everything  together,  (39f)  becomes 

Ri  =Ri  +  7^-urA ,  +  %  ( Vta  ~  vs )  (49a) 

Vhor  K 

=R\  +  DV<a-  dtD  (V,V/;o,  -  RurAs)  (49b) 

R  KVhor 

«Ri  +  ^V,„  -  44  (V,2  -  RrA.)  (49c) 

/y2  \ 

(49d) 

«A-|  4)7  ■)).,!■  (490 


where  we  make  use  of  V2  =  y2  —  RrAs  (13a)  and  14,  ~  V)2  /Vs.  The  latter  is  the  velocity  of 
the  beam  footprint  that  moves  along  the  surface  of  the  earth  and  the  approximation  sign  is 
mainly  due  to  the  fact  that  the  satellite  orbit  is  only  approximately  circular.  Therefore,  the 
first-order  coefficient  of  the  Taylor  expansion  of  Ri{t)  evaluated  at  time  to  can  be  written 
as 


(50a) 

=  (V,r-V,r)-AA-(Vs-Vla)  . 

(50b) 
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Similarly,  we  derive  the  second-order  coefficient  of  the  Taylor  series  expansion  of  Riit) 
by  taking  the  time  derivative  of  (49e),  which  simply  yields  AS  ~  A1] .  Therefore,  the  Taylor 
expansion  of  Ri(t)  (up  to  the  second  order)  evaluated  at  arbitrary  time  to  can  be  written  as 


/?2(0=^o)-^[^o)-'P]  + 


(Vtr  -  Vsr)  ~JUo){Vs~  Vta) 


( t  -  to) 


+  2 


1  (V}~  2VSV, 


svta 


R(to) 


(51) 


Atr  j  (t  —  to)  . 


We  are  now  ready  to  generalize  the  moving  target  range  equation  (5 1 )  for  a  multi-channel 
SAR  system  (i.e.  with  multiple  phase  centers): 


Rp(t)  =  Rq  +  (p  —  l)D  (*¥  —  d>0)  +  (Vtr  -  Vsr)  (t  - 10) 


(52) 


where  Rq  =  R(to),  d>o  =  d>(to);  Vtr,  Vta,  and  Atr  are  defined  for  time  to;  p  —  1,2, 3, ...  (for 
phase  center  1,  2,  3,  etc.);  Vsr  depends  on  Vs  and  d>o  in  a  predictable  way;  Vs,  Vg,  and  Ve 
vary  slowly  with  time  and,  therefore,  may  be  evaluated  anywhere  in  the  neighborhood  of 
to- 


If  we  choose  to  to  be  the  broadside  time  tb,  then  <t>  and  Vsr  ( =  —  Vs<£>)  become  zero,  resulting 
in 


RJt)  =  Rb  +  (p  -  l)DW  +  Vtr  ( t  -  tb)  + 


(P~1)D 

Rh 


C Vta-Vg )  (t-tt) 


1 

+  2 


V} 


2VsVta  .  \  /  ,2 

+  Atr  )  (t  —  tb) 


Rh 


(53) 


where  subscript  b  denotes  the  broadside  time.  We  can  express  (52)  in  a  form  that  is  appli¬ 
cable  to  both  spaceborne  and  airborne  geometries  as  follows: 


Rp(t)  =R0  +  (p-l)D('¥-  d>0)  +  (Vtr  +  Vs<t>0)  Tpri(m  -  m0)  +  ... 


(P~1)D 

Ro 


(vta  -  Vg)  Tpri(m  -  mo)  + 


pn\ 

VsVg  —  2VsVfa  **ir  i  rri  /  w  ,c„  . 

+  —  )  Tpri(m-mo)  (54a) 


2Ro 


Atr 

~2 


Vt 


=R0  +  (p-  1  )D(X¥  -  d>0)  +  VsTpri  2T  _  <$>o  )  (m  -  mQ)  + ... 


(p-VjDVJ) 


pn 


Ro 


Vs 


Vta  Vg  . 

vrvP{m~m)+- 


V2T 2 

s  M  pri 

2Ro 


Vg  ~Vta  R0a  W 

— — 2 - — =Atr  m  —  mo 

vs  vs  v2  IK 


=Ro+D('¥-&o)(p-  1)  +Ds(p  +  ®o)(m-mo)  +  ... 

nn  n  2 

(v  -  p2)  (p  -  1) (m  -  mo)  +  — f  (p2  -  2v  +  rj)  (/w  -  mQy 


Ro 


2R0 


(54b) 


(54c) 
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where 


^  T  7?  \7  r\  — 


‘K2  K 


D, 


Dg  —VgTpri 
Ds  =VsTpri 
Vta 
Vs 
Vtr 
Vx 


V  = 


b  = 


_R0  . 

^  -y2Atr 


m0  =- 


r  S 
*0 


L  pn 


(55a) 

(55b) 

(55c) 

(55d) 

(55e) 

(55f) 

(55g) 


Tpri  is  the  pulse  repetition  interval,  and  m  is  pulse  index  in  slow  time.  Notably,  the  effect  of 
gravitational  acceleration  is  completely  absorbed  into  the  parameter  p2  <  1 .  Clearly,  (54c) 
is  also  applicable  to  airborne  imaging  geometry  by  setting  p2  =  1  since,  for  the  airborne 
case,  the  platform  is  assumed  to  be  moving  along  a  straight  line  and  transmitting  uniformly 
spaced  pulses.  This  assumption,  however,  requires  good  motion  compensation  and  good 
control  of  the  pulse  repetition  frequency  (PRF)  as  a  function  of  ground  speed,  such  that  the 
SAR  ground  speed  Vg  is  kept  the  same  as  its  platform  speed  Vs  (or  Dg  —  Ds). 


Finally,  (54c)  can  be  re-expressed  for  the  case  when  to  =  tj,  as  follows: 

Rp(t)  =Rb  +  Dx¥(p-l)+Ds/j(m-mh)  +  I^l(v-^1)(p-l)(m~mb) 

Kh 

D2  (56) 

+  ^(P2-2v  +  n)(m-^)2  . 

2Kb 

As  stated  earlier,  can  be  compensated  by  channel  balancing  and,  therefore,  may  be 
ignored.  Often,  T]  (the  normalized  target  radial  acceleration)  is  also  ignored  under  the 
assumption  that  the  moving  target  is  not  accelerating.  In  this  case,  (56)  yields 

Rp(t )  =  Rh  +  Dsp(m  -  mb)  +  — —  (v  -  P ~)(P~  1 )  (m  -  mb )  +  x^(P2  -  2v)(m  -  mb)2  , 

Kb  2Kb 

(57) 

where  subscript  b,  again,  denotes  that  the  Taylor  series  expansion  is  evaluated  at  the  broad¬ 
side  time  tb. 


The  accuracy  of  the  equations  of  motion  derived  above  has  been  tested  and  validated  using 
the  recently  acquired  RADARSAT-2  MODEX  data  [5]  [6].  As  discussed,  these  equations 
are  applicable  to  both  airborne  and  spaceborne  imaging  geometries  and  serve  as  a  physical 
basis  for  further  algorithm  development. 
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4  Conclusions 


Both  high  resolution  Synthetic  Aperture  Radar  (SAR)  processing  and  SAR  Ground  Moving 
Target  Indication  (SAR-GMTI)  require  that  a  highly  accurate  imaging  geometry  model 
be  first  established.  This  can  be  quite  easily  accomplished  for  the  case  of  an  airborne 
platform,  which  is  assumed  to  be  moving  along  a  straight  line  and  transmitting  uniformly 
spaced  pulses.  For  an  orbiting  platform,  the  gravitational  force  plays  a  key  role  in  defining 
the  satellite  trajectory  and  its  ground  speed  and,  thus,  must  be  taken  into  account.  Also, 
equations  of  motion  of  a  moving  target  that  accurately  model  a  SAR  system  equipped  with 
multiple  apertures  (e.g.  RADARSAT-2,  TerraSAR-X)  are  evidently  absent  in  the  open 
literature,  partly  because  there  were  no  multi-aperture  spaceborne  SARs  in  the  unclassified 
world  prior  to  late  2007. 

This  short  technical  memo  derives  a  set  of  equations  of  motion  that  accurately  describes 
a  ground  moving  target  in  spaceborne  multi-channel  SAR  imaging  geometry.  The  target 
range  history  model  is  derived  in  the  Earth  Centered,  Earth-Fixed  (ECEF)  frame  of  refer¬ 
ence  using  linearization  for  small  angles  as  a  function  of  the  receive  phase  center. 

The  accuracy  of  the  equations  of  motion  derived  in  this  short  note  have  been  tested  and  vali¬ 
dated  using  the  recently  acquired  RADARSAT-2  MODEX  data.  These  equations  of  motion 
have  been  shown  to  be  applicable  to  both  airborne  and  spaceborne  imaging  geometries  and 
serve  as  a  physical  basis  for  further  algorithm  development. 
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Annex  A:  Acceleration  in  the  ECEF  Frame 


In  this  annex,  we  would  like  to  derive  and  express  the  acceleration  vector  As  of  an  orbital 
satellite  in  the  Earth-Centered-Earth-Fixed  reference  frame.  Let  coe  be  the  earth’s  angular 
rotation  speed,  ignoring  precession  and  mutation.  The  transformation  matrix  T0)  from  the 
Earth-Centered-Inertial  (ECI)  frame  to  the  ECEF  frame  is  then  given  by 

cos(ov)  —  sin(coet)  0 

rro=  sin(coet)  cos  ((0et)  0  .  (A.l) 

0  °  1  _ 

We  begin  our  derivation  with  the  satellite’s  position  vectors  Rs  and  Rv  in  the  ECI  (denoted 
by  the  bar  notation)  and  ECEF  frames,  respectively.  The  transformation  of  the  satellite 
position  vector  between  the  two  reference  frames  can  then  be  accomplished  through  the 
following  two  equations: 

Rj  —  rwRs  (A. 2a) 

R,  =  T^R,  .  (A. 2b) 

Differentiating  (A.2a)  and  (A.2b)  with  respect  to  time,  we  obtain 

v,  =  ravs  +  f  ror^R, 

"0  -1  o' 

=  rroV5  +  coe  1  0  0  R*  (A. 3a) 

_  0  0  0  _ 

v,  =  rJVj+rjR* 

=  +  f  ^rroR, 

"o  10' 

=  r  ^ V 5 + coe  —  1  0  0  R^ 

0  0  0  _ 

Rsy 

—  r^V^  +  cOg  —RSx  ■  (A. 3b) 

0 
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To  get  the  satellite  acceleration  in  the  ECEF  frame,  we  differentiate  again  (A. 3a)  with 
respect  to  time: 


A,  =  r£A,+2f£v,  +  fpl, 

=  r^A, + 2r^r(0v* + 2r£rMR,  +  f^rroR, 


0 

1 

0  ' 

■  1 

0 

0  ' 

-1 

0 

0 

vs  +  (02e 

0 

1 

0 

0 

0 

0 

0 

0 

0 

Vsy  " 

Rsx  1 

T^AiS  +  2coe 

-Vsx 

0 

1 

Q>  ° 

(A.4a) 

(A.4b) 


(A.4c) 


As  is  the  earth's  gravitational  field  (or  acceleration),  which  can  be  approximated  by 


A, 


GMe 

R3s 


R( 


(A. 5) 


where  Me  is  the  mass  of  the  earth  and  G  is  the  gravitational  constant.  This  approximate 
result  may  be  insufficient.  There  are  many  available  models  that  include  higher  spherical 
harmonics  [7]. 
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